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We show that Haldanes new definition of statistics, when generalised to 
infinite dimensional Hilbert spaces, is equal to the high temperature limit 
^— ^ ' of the second virial coefficient. We thus show that this exclusion statistics 

o 

■r-l- ' parameter, g , of anyons is non-trivial and is completely determined by its 

Q\' 

exchange statistics parameter a. We also compute g for quasiparticles in the 

C ' Luttinger model and show that it is equal to a. 
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There has been much interest in the physics of anyons in the past few years. Anyons 
are particles whose many particle wavefunctions pickup a phase e*"'^ under the exchange of 
the positions of any two particles [Q. Arbitary values of a are allowed by the configuration 
space topology of many particle systems in one and two spatial dimensions. In particular 
a = corresponds to completely symmetric wavefunctions (bosons) and a = 1 corresponds 
to antisymmetric wave functions (fermions). The parameter a is traditionally called the 
statistics of the particle, since for a = 1, the antisymmetry of the wavefunction implies the 
Pauli exclusion principle. This has non-trivial consequences for the counting of states in the 
many particle system and hence on the statistical mechanics. 

The effect of the exchange phase on the exclusion principle for arbitrary a was largely 
unexplored until recently. In a seminal paper [0, Haldane has proposed an alternate defi- 
nition of statistics based on a generalised exclusion principle. His definition is for particles 
with finite dimensional Hilbert spaces but can, in principle, lead to the existence of fractional 
statistics in arbitrary spatial dimensions. The parameter governing fractional statistics in 
this case is defined hj g = — ^^^^~ ■, where N is the number of particles and d^ is the 
dimension of the one particle subspace in N - particle Hilbert space. By definition, there- 
fore, (7 = corresponds to bosons and (7 = 1 to fermions. Haldane applied this definition to 
two cases: Quasi-particles in fractional quantum Hall effect (FQHE) systems and spinous in 
quantum antiferromagnets. In both cases he argued that the exclusion statistics parameter, 
g is equal to the exchange statistics parameter a. Johnson and Canright @] have tested 
this prediction numerically for FQHE systems and find this to be true with the modification 
that for quasiparticles, g = 2 — \^ where as a = — |. For the case of anyon gas, Haldane 
argued that since the coupling of particles to the Chern-Simons gauge field cannot change 
the dimension of the Hilbert space, the exclusion statistics of anyons is that of the original 
hard-core particles, ie., (7 = 1. 

In this paper we first generalize Haldane's definition to the case where the Hilbert space 
of the particles is infinite dimensional. We then show that in general the exclusion statistics 
parameter g is proportional to the dimensionless second virial coefficient provided the system 



admits a virial expansion of the equation of state in the high temperature hmit. We apply 

this resuh to the case of anyon gas and find that g = 2a ~ a^, when the exchange statistics 

parameter a , is chosen to be in the range < a < 2. We then analyse the case of anyons in 

a magnetic field, confined to the first Landau level. We obtain g = a for this case. We also 

apply our result to the case of quasiparticles in a one component Luttinger liquid. These 

quasiparticles also have a nontrivial exchange phase due to the braiding properties of the 

vertex operators that create them. Here also we find that g = a, the spinous are a special 

case of this when a = 1/2. Note for quasiparticles with exchange statistics —1/3, we have 

to take a = 2 — 1/3. Thus our definition and calculation is consistent with the numerical 

findings 0. 

The dimension of the Hilbert space of N identical particles when g is the statistics 

parameter is given by, 0,0]. 

(d+{l-g)iN-l))\ 
^^^^^ - Nlid - 1 - giN - 1))V ^^^ 

where d is the dimension of the single particle space. It is easy to see that g = 0(1) yields the 

dimension of the bosonic(fermionic) Hilbert space. It is now straight forward to extract the 

statistics parameter g from eq.(l) in the limit where the single particle dimension d —>■ cxo, 

we obtain 
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g = lim 



m^ - 1 
d^ 



(2) 



2 ^ d^oo N{N - 1 

A regulated definition of the dimension of the Hilbert space is given by the corresponding 
N-particle partition function since we have, D^ = lim^^o ^a^ = liTo^i3-^oTr{e~^^'^), where 
(3 is the inverse temperature and H^ denotes the N-particle Hamiltonian. Therefore, for 
dealing with infinite dimensional Hilbert spaces, we propose a generalisation of Haldane's 
definition. 
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where C is an overall constant of proportionality which we fix next. The definition of g given 
by eq. (3) forms the basis for the rest of the calculations presented in this paper. 



To test the utility of eq.(3), we first consider the case of bosons and fermions in dg space 
dimensions. For convenience we confine all the N-particles in an oscillator potential of the 
form l^(ri,r2, ...,rAr) = \rnuj'^Y^'^i'f'i) where uj is the oscillator frequency and r, are single 
particle coordinates. The potential here merely acts as a regulator and oscillator potential 
which will be used throughout this paper is only a convenient choice. The N-particle partition 
function is then given by the expansion [Q, 
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where the ± signs refer to bosons and fermions respectively. Notice that the second term in 
eq.(4) involves Zi{2j3) since the contribution to this term comes when two particles are in 
the same energy level. Equation (4) is an obvious generalisation of eq.(l) with a cutoff on 
the single particle energies. The single particle partition function is given by, 

-l3Luds/2 

Substituting eq.(5) and eq.(4) in eq.(3) and taking the high temperature limit we find, 
I — (7 = i^T^yaj which essentially isolates the contribution of the subleading order term in 
eq.(4). We therefore set C = 2'^" and immediately obtain g = 0(1) for bosons (fermions). 
Thus factor C occurs because the regulated definition cuts off the energy at 2/? when two 
particles occupy the same state, however for counting purposes we want to cutoff all states 
at the same energy irrespective of their occupation. 

We are now in a position to apply the definition of g as generalised to infinite dimensional 
Hilbert spaces by eq.(3) to N-particle systems with interactions. To begin with consider a 
system of N-interacting bosons in two space dimensions. We first show the connection 
between eq.(3) and the second virial coefficient in the equation of state of the system. The 
virial expansion for the pressure P of a gas is by definition a high temperature or low density 
expansion in terms of the particle number density p = N/V, where N is the particle number 
and V is the volume (area in two space dimensions). In the thermodynamic limit N ^ oo 
and V^ -^ oo when p is held fixed P,^. The virial expansion is then given by 



(3P = p 



1=1 
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where A = J2T^f3/m is the thermal wave length and the dimensionless Bis are called virial 
coefficients which can be expressed in terms of the partition functions. For example the 
second virial coefficient B2 is given by, 

Z2 



B2 = Zi 
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(7) 



(Note: Conventionally B2 has V/X^ as an overall factor instead of Zi when box normalisation 
is used.) Comparing eq.(7) and eq.(3) we immediately obtain g for the case N=2, 

1 
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This is an exact result for N=2. Now consider A^ > 3. The higher virial coefficients are 
given by, 



B3 = ABl - 263, 54 = 952^3 - 165^ + 364. 



(9) 



In general we have, for A; > 3 



Bk = Fk{B2, Bs, ..., Bk-i) + (-1)"-^^ - 1)^. 



(10) 



where Fk are some functions of i?2, -B3, ..., Bk-i whose explicit form is not needed here and 
bk is given by 

[Zi/Zl)^^ 



{rrii} i i 



mA 



(11) 



The summation over m, is constrained by J2^ irrii = k. Thus the virial expansion of the 
equation of state exists iff all b^s are finite. In order to proceed with the computation of 
g for A^ > 3 , let us assume that for the system under consideration the virial expansion 
exists and hence all -B^'s are finite. The high temperature expansion for the factor Z^/Z^ 
in eq.(3) is in general given by 
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where /2 is some function of the interaction strength and may on N. In the above equation 
we have only shown the expansion up to next to leading order as /3 — > 0; higher order terms 
are irrelevent for our purposes. Using eq.(12) we immediately see that the leading divergence 
in 6fc as /? — s> comes from the overall factor Z^"^ and is given by \ j [^j5uS)'^'^^~^^ ■ For 6^ to 
be finite in the limit /? — i> 0, we must demand that all terms in the sum involving powers 
of (Zj/Z^) in eq.(ll) up to order {^(5lJJ)'^-'^^^^^^^^^ must cancel order by order. This necessarily 
implies that the coefficient of {P^Y' ^lust be set to zero. Demanding this we immediately 
find, 

Af-2 AN-n) A2) 

/f ) = Y. ^i-l^-'^ = N{N - 1)1^, (13) 

n=l ^- -'^ • 

where all /2 are now related to /2 corresponding to the two particle partition function. 
Substituting eq.(12) and eq.(13) in eq.(3) we obtain for g, 

\-9 = f?^ = -2i?2 (14) 

We have, therefore, now a relation between g and the second virial coefficient B2 for a 
system for which the virial expansion exists. Notice that in obtaining the result (14) for 
g, it is sufficient to demand that order (/Jcu)^ terms cancel in the sum given in eq.(ll). 
This is obviously a much weaker condition than demanding that the virial coefficients be 
finite. Nevertheless the whole analysis makes sense only if the virial expansion is valid. 
This then is the main result of the paper: The Haldane exclusion statistics parameter g 
is completely determined by the high temperature limit of the second virial coefficient of a 
system of interacting particles which admits a virial expansion in this limit. 

We can now apply this result to some well known systems. We first consider anyon gas. 
Here the second virial coefficient is well known and is given by [^, 

B2 = ^ , (15) 

where a = is bosonic and a = 1 is fermionic. Using eq.(14) we therefore have, 

g = a{2-a) (16) 



which has the correct hmit for a = 0, 1. This result would be true for two anyons. However, it 
is not conclusively proved that the higher virial coefficients are finite in the case of anyon gas. 
It has been proved that the third virial coefficient is finite @J^. If indeed all virial coefficients 
are finite, then the result given in (16) would be true for anyon gas in general. This is an 
interesting result since the exclusion statistics parameter g is not the same as the exchange 
statistics a unlike the systems considered by Haldane. It was expected that the particles 
constituting anyon gas would be classified as hardcore bosons (hence fermions) since the 
coupling of a particle to Chern-Simons gauge field does not affect Hilbert space dimensions. 
Obviously our result does not bear out this argument. It is easy to see why. The two anyon 
spectrum, after taking out the centre of mass, is given by E^^i = uj{2n + |/ — a| + 1). When 
a is nonzero all the energy levels are shifted up by a (/ < 0) or 2 — a; (/ > 0) ( equivalently 
—a) no matter where the cutoff lies. 

Next consider anyons confined in a magnetic field. Once again g is given by eq.(16) 
when the partition function includes the trace over all Landau levels. This is easy to see 
since both the oscillator frequency u and the cyclotron frequency Uc act as regulators and 
at high temperatures one recovers the anyon gas result. If however we confine the trace to a 
single Landau level, the we obtain g = a in conformity with the results obtained for FQHE 
systems. 

We will now consider quasiparticles in the Luttinger model. The low energy physics of 
this model maps on to the massless Thirring model. As is well known 0, the model is 
exactly solvable. We will therefore be able to compute the RHS of equation(3) directly in 
this case. The theory can be written completely in terms of the left(L) and right (R) currents 
satisfying the algebra [J^, J^] = nSmn^rs, where n,m = 1,2,... and r,s = L,R. We also 
have the zero modes [0q, JJ] = i6rs- The hamiltonian is 
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Quasiparticles are created by the action of the vertex operators V^{x) =: e *9'^<^'^(^) : on the 
ground state. The 0'"(x) are the bosonic phase fields given by. 
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where a^^R = ±1. The periodicity properties of the compact zero modes 0o' constrain the 
allowed values of qr to be qr = {NR + o"r^) and R is the radius of the field 0. It is a 
function of the interaction strength. Our conventions correspond to i? = 1 for the case of 
non interacting fermions. The constraints on g^ imply that quasiparticles in the left sector 
must be created along with quasiparticles in the right sector (except for special values of 
R). However because the hamiltonians of the two sectors are completely decoupled, it is 
consistent to analyse the spectrum of the two sectors independently. We will therefore focus 
on the left sector alone. 

We first consider the space of one quasiparticle states which we define to be the space 
of the states \x >= Vl{x)\0 >. These are not a linearly independent set. From the form of 
Vl{x), it follows that \x + L >= e^°'^\x >, where a = q\. Therefore, we have the expansion 
\x >= X]^ie*^"''|n >, where A;„ = ^(n + |). From the form of Vl{x)\0 > it follows 
that n > 0. It can be easily shown that \n > form an orthogonal set of eigenstates of the 
Hamiltonian with eigenvalues E^^ = vpkn- The single particle partition function is then 

given by, 

p-/3"/2 - , ~ 27r 

Z, = ^ ~ = e^^'^I^Zf- a = '-fpv,. (19) 

1 — e~^ -tv 

Next we come to the N quasiparticle states which we define to be the span of \{xn} >= 
n^=i1/2(x„)|0 >. From the fact that Vl{x)Vl{y) = e'''''Vl{y)Vl{x) , it follows that these 
many quasiparticle states pick up a phase e*"'^ under the exchange of two of the coordi- 
nates. We normal order the vertex operators to obtain, \{xn} >= I\n>m 
'[\n=i ^M^n) • |0 >• Again from the form of the vertex operators it follows that we have the 
expansion, 

E'^^({UIK})|{U>, (20) 
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where (t)%{{kn}\{xn})i are symmetrized N particle plane waves with momenta {kn}-, where 
kn = ^{iT- + ^)j n > 0. The states \{kn} > can be shown to be eigenstates of the 



hamiltonian with energy E^{{kn}) = VFJ2n=i ^n- We can show that |{A;„} > are a hnearly 
independent set of states. The states with the same energy are however not orthogonal. 
The N-quasiparticle partition function is then Zjy = e~^^ °'^'^Z^, where Z§ is the N-particle 
bosonic partition function. We can now exactly compute g using the high temperature ex- 
pansion for Z^ and obtain g = a. Thus the exchange and exclusion statistics parameters are 
identical for these models. Note that the exchange statistics in one dimensions is somewhat 
arbitary. We could have changed it by multiplying the vertex operators with suitable cocy- 
cle factors. The exclusion statistics parameter is however unambigous and unique. When 
R = -y=, the theory is equivalent to the low energy physics of the SU(2) symmetric quan- 
tum antiferromagnetic chain. We then have the spinon excitations with a = 1/2. Thus we 
recover Haldane's result that g = 1/2 for this case. 

The above example gives a clear insight into the mechanism of the phenomenon. What 
is happening is that the addition of a quasiparticle causes a phase shift of every other 
quasiparticle, resulting in an energy shift of ^^j^ per particle . When we count the dimension 
of the single particle space with a fixed (smooth) cutoff, there are a states missing. The 
important thing here is that the all single particle levels shift up by the same amount however 
high the energy. This is why we get g to be well defined and nontrivial in the cutoff going 
to infinity limit. 

Exactly the same thing happens in the case of the anyon gas except that the relevent 
levels do not all shift by the same amount. The / = 0, —2, —4, ... levels in the two anyon 
spectrum shift up whereas the positive / = 2,4,6,... levels shift down by an amount —a. 
Equivalently the later set may be considered as made up of levels / = 0,2,4,6, ... with an 
upward shift given by 2 — a. This is why g is not equal to a but has a nonlinear dependence 
on it given by the product of the two energy shifts, namely a{2 — a) as though the full g 
is the product of g^s corresponding to two one component systems. However the important 
fact is that g is completely determined by a. When we consider anyons in a strong magnetic 
field and restrict ourselves to the first Landau level, then again we have Vs of only one sign. 



This again results in g = a. 

Thus the fundamental character of exclusion statistical interactions seems to be that 
they cause scale-invariant phase shifts and hence scale-invariant energy shifts. How this is 
connected to the non-orthogonality of position eigenstates which is stressed in ref. is not 
clear to us. We also note that these scale-invariant phase shifts, if they occur in any model, 
is a non-perturbative effect and the model cannot be analysed by perturbation theory JlOl . 



Thus a non-trivial g implies the inapplicability of standard many-body perturbation theory, 
as has been stated in ref 0. 

Finally, we have related g to the high temperature limit of a thermodynamic quantity, 
which in principle, is much easier to measure than the exchange phase of quasiparticles in 
condensed matter systems. It may therefore be possible to use this connection to device a 
realistic experiment to measure the Haldane exclusion statistics of quasiparticles. 

We gratefully acknowledge many interesting discussions with G.Baskaran. 
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